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We present numerical results for the dissociation cross 
radially-excited charmonia in collisions with light mesons, 
ativistic quark model, so all parameters are determined by 
Examples of dissociation into both exclusive and inclusive 
tion cross sections of several C=(-|-) charmonia may be of 
heavy ion collisions, since these states are expected to be 
The relative importance of the productions of ground-state 
a pion-charmonium collision is demonstrated through the 
cross sections. 



sections of ground-state, orbitally- and 
Our results are derived using the nonrel- 
fits to the experimental meson spectrum, 
final states are considered. The dissocia- 
considerable importance for the study of 

produced more copiously than the J/f/). 

and orbitally-excited charmed mesons in 

V^-dependent charmonium dissociation 



I. INTRODUCTION 

The interactions of mesons in systems containing both 
light and heavy quarks have long been of interest to 
hadron physicists. For example, some models predict 
that the open-charm D and D mesons have sufficiently 
attractive residual strong interactions with nucleons to 
form "charmed nuclei" 0. Charmed hadron bound 
states may exist in other systems as well; Novikov et 
al. long ago speculated that the nominal 3^ Si cc state 
i/'(4040) might actually be a quasinuclear "molecule" 
bound state of a D* D* pair Several quark model 
studies have shown that Q^q^ mesons should exist for 
sufHciently large heavy-quark mass toq (for a recent re- 
view see Richard *^). In a recent "pedagogical" applica- 
tion, the toq — > oo limit of the Q^if' heavy- light system 
(the so-called BB system) has been used as a theoretical 
laboratory for the study of nuclear forces, and nuclear po- 
tential energy curves have been derived using the nonrel- 
ativistic quark model and lattice gauge theory (LGT) 
i 

Recently, further interest in the interactions of light- 
and heavy-quark mesons has arisen in the context of 
heavy-ion collisions and the search for the quark-gluon 
plasma (QGP). One signature proposed for the identi- 
fication of a QGP is the suppression of the rate of 
formation of the J /tp and other cc bound states. The 
long-ranged linear confining potential between a cc pair 
would purportedly be screened by a QGP, so a cc pair 
produced in the collision would be more likely to separate 
than to populate bound cc resonances. 

Direct experimental confirmation of such a suppression 
can be detected for example through the observation of 
lepton pairs from the decay J/V' — ^ The simplest 

interpretation of an observed J /ip I'^t' signal would 
be to assume that all J/V' mesons survive until they decay 
outside the interaction region. However, if dissociation 
reactions such as tt -I- 3 jip — > D*D and p + J/^ DI) 



are important, the interpretation of the experiment is 
more complicated; a weak J jil) signal might simply be 
due to dissociation through such "comover absorption" 
processes. The actual size of these low-energy charmo- 
nium dissociation cross sections is currently very contro- 
versial, and their evaluation is the subject of this paper. 



II. APPROACHES 

A very wide range of theoretical estimates of low- 
energy charmonium dissociation cross sections has been 
reported in the literature, largely due to different as- 
sumptions for the dominant scattering mechanisms. We 
will briefly review the three main approaches used be- 
fore presenting new results from our quark model calcu- 
lations. 



A. Quark Interchange 

Quark model calculations of charmonium - light 
hadron cross sections were first reported by Martins, 
Blaschke and Quack |8|, who used a quark-interchange 
model 9] to treat tt-I- J/-0 collisions. This reference used 
standard quark model one gluon exchange (OGE) forces 
(spin-spin hyperfine and color Coulomb), augmented by 
a color-independent confining force that was assumed to 
act only between qq pairs (no qq or qq anti-confining in- 
teraction). Since the color Coulomb terms experience 
destructive interference between diagrams (due to color 
factors) and the large mass of the charm quark makes 
the hyperfine term rather weak, Martins tt al. con- 
cluded that OGE forces alone give rather small cross 
sections; they estimated a cross section for tt -|- J/ip 
D*D + h.c. + D*D* of « 0.3 mb at « 4.2 GeV. 
However their color-independent, gg-only confining in- 
teraction had no color factor cancellation, and so gave a 
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much larger peak cross section of « 7 mb near 4.1 GeV. 
These two exclusive processes {D*D + h.c. and D*D*) 
were found to peak quite close to threshold, and then fell 
rapidly with increasing invariant mass due to suppres- 
sion from the Gaussian meson wavefunctions assumed. 
(See Fig.2 of Rcf. 8].) Subsequently, Wong et al. [13 
applied the same approach to this problem, albeit with 
the conventional quark model A • A dependence for all 
terms in the interquark Hamiltonian, including the lin- 
ear confining potential. The interaction terms assumed 
were again color Coulomb, spin-spin hyperfine and lin- 
ear scalar confinement. The destructive interference be- 
tween all diagrams, due to zero-sum color factors, led to 
a rather smaller tt -|- J/^p —* D*D + h.c. + D*D* cross 
section, with a rather broad peak of roughly 1 mb near 
4.1 GeV. Wong et al. also considered ir + ip' , p + J/4> and 
p + ip' dissociation, and found that these cross sections 
are much larger than n+J/ip near threshold, due to more 
favorable kinematics. The p-\- (cc) processes are exother- 
mic, and so actually diverge at threshold. Finally, Wong 
et al. studied the importance of the so-called "post-prior 
ambiguity" in these calculations; the use of exact 
qq Hamiltonian wavefunctions in the present paper elim- 
inates much of this systematic effect, but an important 
discrepancy remains due to the use of relativistic phase 
space and physical masses. 

Shuryak and Teaney .13j gave a comparable rough es- 
timate of « 1.2 mb for a low-energy TT + J/tp cross section 
driven by the nonrelativistic quark model's spin-spin in- 
teraction. Actually the specific process they considered, 
TT+J/ij] — > Vc + P, is zero at Born order due to a vanishing 
color factor; color was not incorporated in their estimate. 



B. Meson Exchange 



was much larger than the D-exchange results of Matinyan 
and Miiller, due to new three- and four-meson vertices in 
their effective lagrangian. Haglin introduced a sim- 
ilar SU(4) symmetric meson lagrangian, and also found 
rather large cross sections of 5-10 mb for y/s — 4-6 GeV 
for many low energy charmonium dissociation reactions 
(see for example Fig.2 of Ref . [Til ) . Subsequent work by 
Haglin and Gale 19] showed that the tt+J/4> total inelas- 
tic cross section would reach an extremely large value of 
roughly 100 mb at ^/s = 5 GeV, and p + J/ip a fantastic 
« 300 mb (with both still increasing) in this model, as- 
suming pointlike hadron vertices. Similar large tt -|- J/ijj 
and p + J/ip cross sections have been reported by Oh, 
Song and Lee in pointlike meson exchange models [23| . 

Navarra et al. f2l| have recently questioned the as- 
sumption of flavor SU(4) symmetry; keeping only isospin 
symmetry, they find rather smaller cross sections for 
TT + J/^ D*D + h.c, ca. 20-25 mb at = 5 GeV. 
They confirm that D* exchange is much more important 
numerically in meson exchange models than the D ex- 
change originally assumed by Matinyan and Miiller. 

Of course it is also incorrect to assume pointlike hadron 
form factors. This has been noted both by Lin and Ko 
[l7l | and Haglin and Gale Both collaborations in- 

vestigated the effect of assuming dipole forms for the ef- 
fective three-meson vertices, and concluded that the pre- 
dicted cross sections were greatly reduced (once again to 
typically 1-10 mb scales) with plausible vertex functions. 
(See for example Fig.4 of Ref.flT^ and Fig.7 of Ref.[l|.) 
Accurate calculations of hadronic vertex functions are 
clearly of crucial importance for meson exchange mod- 
els of charmonium dissociation. Some results for these 
form factors, obtained from QCD sum rules, have been 
published by Navarra et al. |^ 0| . 



Charmonium - light hadron scattering can also take 
place through t-channel meson exchange. This mecha- 
nism was first discussed by Matinyan and Miiller [l^ . 
who were motivated to study tt -t- J inelastic scat- 
tering by the great discrepancy between the ca. 7 mb 
quark-model result of Martins et al. and the very 
small low-energy cross sections found using the Peskin- 
Bhanot approach [l5l |. In the meson-exchange picture, 
charmonium dissociation reactions proceed through t- 
channel exchange of charmed mesons such as D and D* . 
Matinyan and Miiller assumed only D exchange, and 
found mb-scale cross sections for the two low-energy dis- 
sociation processes TT -I- J/-)/) — > D*D + h.c. and p + J/4> — ^ 
D*D*. 

This work has since been generalized to other t-channel 
exchanges and effective meson lagrangians. Lin, Ko and 
Zhang had previously proposed an SU(4) flavor symmet- 
ric vector-pseudoscalar meson effective lagrangian which 
they had applied to open-charm meson scattering [T^. 
Application of this same lagrangian to the ir + J/tp disso- 
ciation reaction ir+J/ip — > D*D + h.c. gave a rather large 
cross section of « 20-30 mb for ^ = 4-5 GeV 0; this 



C. DifFractive Model 

A high-energy diffractive description of scattering of 
heavy quarkonia which was developed in 1979 by Peskin 
and Bhanot li?] has also been applied to the calcula- 
tion of charmonium cross sections. It should be stressed 
that this method is only justified at high energies [25| . 
and then only for deeply-bound QQ systems. It is in 
essence a gluon-sea model of high-energy diffractive scat- 
tering of physically small, high-mass Coulombic bound 
states by light hadrons. This model predicts reasonable 
mb-scale cross sections for J/ijj hadronic cross sections 
at > 10 GeV 0, IH I2I i^. At low energies, 
however, this mechanism taken in isolation predicts ex- 
tremely small (sub-/ib) cross sections for J/ip + tt and 
J/^p + N (see Fig.2 of Ref.'2l). Presumably this means 
that the Peskin-Bhanot diffractive scattering mechanism 
is unimportant in the low-energy regime of greatest rel- 
evance to QGP searches, and other mechanisms such 
as quark interchange and meson exchange dominate cc 
strong interactions at these low energies. Indeed, a re- 
cent comparison |28j | with lattice gauge computations 
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shows that the operator product expansion breaks down 
for quark masses below roughly 10 GeV, and therefore 
the Peskin model is inapplicable to light and charmed 
hadronic physics. 

Redlich et al. |29| have argued that these diffractive 
cross sections are actually accurate at low energies, and 
if combined with vector dominance can account for the 
experimental — > + iV — > open charm cross sec- 
tion, whereas the much larger meson-exchange result for 
a{J/ip + N ^ D + Ac) of HagHn greatly exceeds 
the charm photoproduction measured in experiment, as- 
suming vector dominance through the J/V'- Hiifner et 
al. |2g however argue that this test is misleading, as the 
assumption of vector dominance through the J/tp alone 
is unrealistic for these processes. Since the ip' is much 
closer to open charm threshold than the J/V', and is pre- 
dicted to have larger dissociation cross sections to open 
charm, this is clearly a potential source of inaccuracy for 
any J/ip-oidy vector dominance model. 

Finally, QCD sum rule calculations have been per- 
formed 30] which find a much larger low-energy cross sec- 
tion than the diffraction model Ref. [Tsj l , and are in rough 
agreement with quark interchange results near threshold. 
Although there is approximate agreement of scale at low 
energies, the sum rule results find that the exclusive cross 
sections increase monotonically with energy. We believe 
that this is incompatible with hadronic form factors, 
which may require the inclusion of higher-dimensional 
operators in the sum rule calculations. 



D. Synopsis 

Clearly the scale of charmonium dissociation cross sec- 
tions at low energies remains an open question. Nei- 
ther the experimental values nor the dominant scattering 
mechanisms have been convincingly established. In this 
currently rather obscure situation we can best proceed 
by deriving the predictions of the various models and 
searching for the least ambiguous comparisons with ex- 
periment, in as unbiased a manner as possible. Here we 
attempt to contribute to this research through a care- 
ful and detailed study of the predictions of one of the 
theoretical approaches, the quark interchange model. 



III. QUARK INTERCHANGE MODEL 

The Born-order quark interchange model approximates 
hadron-hadron scattering as due to a single interaction 
of the standard quark-model interaction Hamiltonian Hj 
between all constituent pairs in different hadrons 9]. In 
the current study we specialize to quark interactions that 
are simple potentials times spin and color factors. 

Hi = (^vcou. (r) I + Vconft. (r) I + v,s (r) S,; • S,^ T"^ • TV 

(1) 



The potentials are the standard quark model color 
Coulomb, linear confinement, and OGE spin-spin hy- 
perfine terms, vcou. = a^/r, Wcon/t. = -36r/4, v^s = 
— {Swas /dmiUij) Sa-{r). (The Gaussian-regularized delta 
function is Sa{f) = a^/n^^^ ■ e~'' ^ .) 

Since this Hamiltonian is T° • in color space, quark 
line rearrangement is required to give nonzero scattering 
amplitudes between initial and final color-singlet hadrons 
at leading order in Hj. In the case of qq meson- meson 
scattering, this Born-order amplitude is given by the sum 
of the four "quark Born diagrams" shown in Fig.l. Each 
interaction in each diagram has an associated "signature" 
fermion permutation phase, color factor C, spin matrix 
element S, and a spatial overlap integral 2. The eval- 
uation of these various factors is discussed in detail in 
Ref. 9j. There are several simplifications in qq meson- 
meson scattering; the signature phase is always (—1), the 
fiavor factor is diagram independent (and is unity here), 
and the color factors are (—4/9) (capture) and (-1-4/9) 
(transfer). The full meson- meson T-matrix is given by 
the sum of color Coulomb, linear confinement and OGE 
spin-spin T-matrix elements, each of the form 

tAb^cd ^ (_i).jr./(_4/9).(5^j)^^ + (_4/9).(5^j)^^ 



+ (+4/9) • (5 ® T)ti + (+4/9) • (5 ® T)t, I . (2) 

The angle brackets refer to the fact that the spin- 
and space-matrix elements do not always factor, 
and must in general be evaluated using a Clebsch- 
Gordon series. (This complication applies to spin- 
triplet, orbitally-excited mesons.) To evaluate the 
cross section for a given reaction at a given energy, 
we first evaluate the overlap integrals (given below) 
for each set of orbital magnetic quantum numbers, 
{Lc, Lcz\Ld, Ldz\ I \La, Laz\ Lb, Lbz), using an adap- 
tive Monte Carlo technique. In this method we fix 
A — Az, so the overlap integrals are functions of Sic- 
(The magnitudes of A and C are determined from ^/s 
and the physical meson masses using relativistic kine- 
matics.) We then evaluate spherical harmonic moments 
cim = J <^^C^m(^c)^(^c) of the Overlap integrals, for 
each diagram and interaction, usually up to Z = 4. These 
spatial overlap integrals are then combined with the spin 
matrix elements {Sc, Scz] Sd, Sdz\ S \Sa, Saz] Sb, Sbz) 
of I and Si • S^ in a Clebsch-Gordon series to form the full 
T-matrix element 

rpAB^CD _ I^J^^ J^^. J^^ J^^l rp I j^^ j^^. j^^ j^^^ ^ 

Polarized cross sections are then given by 

AB^CD 4:EaEbEcEd I Pel f ,rpAB^CD\2 

= ^ j dnc IT,. I 

(4) 
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FIG. 1: The four quark-interchange meson-meson scattering 
diagrams in the "prior" formalism. 



and the unpolarized cross sections given here are deter- 
mined by summing over magnetic quantum numbers as 
usual. 

Since our Hj consists of simple potentials, it is con- 
venient to evaluate the overlap integrals in real space. 
(In most recent studies we evaluated these overlap in- 
tegrals in momentum space, since they can be conve- 
niently expressed as convolutions of the quark-quark 
momentum-space T-matrix with external meson wave- 
functions; see for example Ref.^3|- The real-space over- 
lap integrals, which may be obtained by introducing 
Fourier transforms in Eqs.(l-4) of Ref.[lQ| for the spe- 
cial case Tft{q,pi,p2 ) ^Tf,{q), are 



Jci = / <fx tpAixA) '4'b{xb) V'c(^) i'*DixA + Xb ~ r) 



v{r) exp<{ -^(A + ncC) ■xa--^{A~hdC) -XB+iA-f 



(5) 



/ 



= / d^X t/jAixA) ^b{xb) i^hi^A +XB ~r) Voi^) 



v{r) exp | + ^(A-^cC') ■xa-^{A- ^idC) -fs+iA-r | 

Where the measure is d^x = (Pr (Pxa'I'^xb /{2n)^. We 
also introduced = 2mq/{mq + rriq), and the identities 
MA = /^B = 1 and /ic + fj-o — 2 were used in deriving the 
overlap integrals; these relations are valid for processes of 
the type {nn) + {cc) — > {nc) + {cn). The spatial wavefunc- 
tions above are the usual nonrelativistic quark potential 
model functions ip{rqq), normalized to J d^r \tp{r)\'^ = 1. 
The wavefunctions employed in this paper to evaluate 
these overlap integrals are numerically determined eigen- 
functions of the full quark model Hamiltonian (with in- 
teraction given by Eq.(l)). 

Since we use relativistic phase space and physical 
masses in evaluating our cross sections, there is a post- 
prior ambiguity in our results • The overlap integrals 
given above are the "prior" form, in which the Hj inter- 
action takes place prior to rearrangement (Fig.l). In the 
"post" form, rearrangement followed by interaction, the 
scattering amplitude is given by a different set of spin 
matrix elements and overlap integrals. The post overlap 
integrals for the two capture diagrams are 



d\ i>A{f) TpBixc +XD -r) ipc{xc) i>*DixD) 



v[r) exp<j +^[A + iiDC)-xc-]^{A-^jLDC)-XD-iC-r 



jpost ^ / jg^, i/,^(f^ _|_ _ f ) il}B{r ) ^Pci^c) i^hixo) 



(9) 



v(r) expl^ + ^{A-pcC)-xc-^{A+ficC)-XD+iC-r^ 

(10) 

The transfer diagrams Tl and T2 (Fig.l) in post and 
prior formalisms are identical. 



IV. RESULTS 



, , / * (a C) " * (a C) 'A ^1 have obtained Born-order quark-model results for 

v(r) exp < + ^ ( fic ) XA+ ^[ +fJ'D ) xb i r> j) dissociation cross sections of ground-state and excited 

(6) charmonia into exclusive final states, and ii) total inelas- 
tic dissociation cross sections from the quark interchange 

/'Q , ,^ \ ^ ,^ \ ,^ mechanism. These results assume the quark interchange 

d X ^a(xa) Vs(xs) ^c{xb +r)i:o{xA~r) ^^^^^ described in the previous section. In the quark 

model Hamiltonian we assume a quark-gluon coupling 
constant Ug = 0.6, a — 0.9 GeV in the spin-spin hypcrfine 
J i I X r'\ ^ ^ ( A n\ ^ ■ A ^\ term and a string tension of 6 — 0.16 in the linear con- 
v{r) el^^^--[A + ^lcC)■XA + ^[A + ^lDC)■XB+^A■r finement term. The light and charmed quark masses are 

(7) taken to be 0.33 GeV and 1.6 GeV, respectively. These 
parameters can reproduce the 1 = 2 tttt S'-wave exper- 

T /"is if^\ir^\i*f^ ^^;*/^ imental phase shifts and are used in the following 

It^= dx tpAixA) i^B[XB) ^cK^B +r)^jj[xA-r) subsections to calculate charmonium dissociation cross 
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sections which are average values obtained with "prior" 
and "post" forms. Cross sections for 7r+ scattering are 
presented in the following sections, other pion cross sec- 
tions may be obtained assuming isospin symmetry. 



A. Dissociation cross sections of excited charmonia 

The principal mechanism for production of charmonia 
at small-x in heavy ion collisions is thought to be the two- 
gluon fusion process gg — > cc. One therefore expects J /if) 
production to be relatively weak, since the formation of 
C=(— ) states requires an additional gluon. The C=(-|-) 
mesons that have especially large gg couplings, such as ric, 
Xco and (to a lesser extent) Xc2 and their, as yet uniden- 
tified, radial excitations should be the dominant cc states 
produced. The relative strengths of cc couplings to glue 
are dramatically illustrated by the total widths of char- 
monia; the r^c and XcO total widths, thought to be due 
mainly to cc — > gg, are two orders of magnitude larger 
than the J total width. Even the smaller Xc2 width is 
roughly 20 times the J/^ width. This suggests that the 
production of charmonia from a quark gluon plasma is 
probably dominated by these C=(-|-) states rather than 
J/V', so the evolution of C=(-|-) states produced in a 
heavy-ion collision may be more important for under- 
standing charm production than the J /tp- The possibil- 
ity that much of the J lij) signal originates from radiative 
transitions of parent XcJ states (Sll |33 also suggests that 
an understanding of the interactions of these C=(-|-) and 
excited cc states with light hadrons may be of great im- 
portance for simulations of heavy flavor production in 
heavy ion collisions. 

It is straightforward to determine the dissociation cross 
sections of cc states other than the J /tp in the quark in- 
terchange model; one simply changes the external state 
attached to each of the four scattering diagrams of 
Fig.l. There is a technical complication with spin-triplet, 
orbitally-excited charmonia, since the spin and space de- 
grees of freedom do not factor trivially in these states, 
unlike the scattering of S-wave mesons considered previ- 
ously jllj . Instead we must evaluate overlap integrals and 
spin matrix elements for each set of magnetic quantum 
numbers, which are then combined using the appropriate 
Clebsch- Gordon coefficients to give scattering amplitudes 
of mesons with definite J (e.g. tt + Xc.j DD*). 



B. Total TT-charmonium dissociation cross sections 
from constituent interchange 

We have evaluated the ttJ/i/j, TTip' , and ttXcJ exclusive 
and total cross sections up to 4.5 GeV in the center of 
mass frame. The six final states with nonzero couplings 
in the model which are open in this regime are 

TT+ ^ cc ^ D° D*+ + cc. (11) 



7t J/\|( total 




3.8 3.9 4 4.1 4.2 4.3 4.4 4.5 



FIG. 2: Theoretical -k^ J/ip cross sections in the quark in- 
terchange model. The figure shows all nonzero partial cross 
sections open to ^/s = 4.5 GeV; the total cross section, ob- 
tained by summing these and their charge conjugate channels, 
is shown as a solid line. 




3.8 3.9 4 4.1 4.2 4.3 4.4 4.5 



FIG. 3: IT XcO cross sections, legend as in Fig. 2. 

7r+ + cc -> (12) 

7r+ -h cc + cc. (13) 

7r+-f cc^L)°Di+(^Pi)-Hc.c. (14) 

■K+ + cc D° D*+ + cc. (15) 

■K+ + cc^ D*^^D+{^Pi)+c.c (16) 
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n\, total 




FIG. 4: TT Xci cross sections, legend as in Fig. 2. 
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FIG. 5: 7r"'"Xc2 cross sections, legend as in Fig. 2. 



The masses of the D*+ , D*° , D+ and D*+ mesons 
are taken from the 2002 PDG compilation |3l|. The 
vector meson is assumed to have a mass of 2.5 GeV. 
The spin-triplet state Dl~^{^Pi) and the spin-singlet state 
D^{^Pi) mix to form the observed state £'i(2420)"'" and 
another unobserved state [s^- We assume masses for 
the Dl+{^Pi) and D+{^Pi) of 2.427 and 2.4 GeV, re- 
spectively. 

The total inelastic cross sections to 4.5 GeV are shown 
in Figs. 2-8. Results for exclusive reactions are as shown 
in the figures. The total cross section also includes charge 
conjugation final states where appropriate. 

We note the following general features of the cross sec- 



ic Y total 




FIG. 6: n'^ip' cross sections, legend as in Fig. 2. 



p J/y total 




3.8 3.9 4 4.1 4.2 4.3 4.4 4.5 



FIG. 7: p V' cross sections, legend as in Fig. 2. 



tions. All cross sections rise rapidly according to thresh- 
old kinematics and subsequently fall off at a scale of 
^QCD as expected for exclusive flavor exchange reactions. 
The ip' cross sections are roughly ten times larger than 
corresponding ip cross sections. This is in accord with 
the ratio r'^Jr^p w 2 and the notion that cross sections in- 
crease with hadron size, although we stress that no simple 
scaling relationship exists in the quark interchange model. 
We note that the ratio of -0' to "0 is substantially smaller 
than the factor of 5000 predicted by the Peskin-Bhanot 
computation [23| ■ The cross sections tend to fall more 
rapidly than those for ip. We suspect that this is due to 
the node in the ■0' radial wavefunction which manifests 



7 




4.4 4.45 4.5 4.55 4.6 4.65 



FIG. 8: p^tp' cross sections, legend as in Fig.2 but continued 
to^/s = 4.65 GeV. 

itself as a zero in the cross section about 200 MeV above 
threshold, causing the cross section to drop more rapidly 
than that for the ground state. Again, this feature can 
be expected to be quite general. Finally, we note that 
the XcJ cross sections grow with angular momentum J, 
which is naively expected due to the increasing number 
of Jz states present. 

A simple parameterization of these cross sections may 
prove useful for further numerical investigations. We 
have found that it is possible to fit many of our numer- 
ical cross sections with a simple functional form which 
is motivated by the expected threshold behavior with an 
exponential decay representing suppression due to flavor 
exchange: 



O-(s) = CTmax ( I exp {p{l - c/emax)) (17) 

where e = ^/s - Mc - Md and p = 1/2 + L^f^ for 
endothcrmic reactions and p = — 1/2 + L*^^ for exother- 
mic reactions. Here Lmin refers to the minimum possible 
value for the total orbital angular momentum of the final 
state consistent with conservation of angular momentum 
and parity. One expects this to dominate the thresh- 
old behavior of a given reaction. In practice we do find 
that many cross sections are well described by assuming 
that the orbital angular momentum in the initial chan- 
nel is zero, however in general many waves contribute 
and it is more convenient to simply fit the value of p. 
We have found that this procedure describes all of our 
■n + cc reactions quite well (however, this is not true for 
p + cc). Results for the parameters p, (Jmaxi and Cmax are 



TABLE I: Cross Section Fit Parameters 
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0.026 


0.127 


^max 


0.059 


0.044 


0.074 


0.050 


0.050 


0.052 


TT*' -» 














P 


0.67 


1.17 


0.84 


0.74 


0.83 


1.42 




13.04 


6.30 


2.66 


1.78 


0.466 


3.81 


^max 


0.027 


0.034 


0.033 


0.032 


0.032 


0.037 


TTXcO -> 














P 


1.60 


1.86 


1.34 


1.59 


0.90 


2.03 


max 


2.02 


0.555 


0.380 


0.244 


0.054 


0.306 


^max 


0.147 


0.114 


0.123 


0.130 


0.141 


0.108 
















P 


1.63 


1.84 


1.44 


0.95 


1.57 


1.93 


(^max 


3.02 


1.03 


0.576 


0.406 


0.092 


0.578 


^max 


0.154 


0.122 


0.130 


0.120 


0.132 


0.114 


1"Xc2 -> 














P 


1.65 


1.84 


1.10 


1.52 


1.36 


1.89 


(^m.ax 


3.64 


1.37 


0.724 


0.598 


0.122 


0.790 


^max 


0.157 


0.127 


0.132 


0.138 


0.133 


0.116 



presented in the Table I. The and Di referred to in 
the column headings represent the Di{^Pi) and £'i(^Pi) 
states respectively. All channels except D*D* include 
charge conjugate reactions in the fit parameters. 



V. SUMMARY AND CONCLUSIONS 

Total charmonium dissociation cross sections have 
been computed up to 4.5 GeV in the center of mass of 
the TT + cc system. The computations employ standard 
constituent quark model dynamics and all parameters 
arc fixed by spectroscopic data. Exact numerical wave- 
functions have been employed in the computations to 
minimize post-prior discrepancy. We have also presented 
results for positive charge conjugation Xc.i dissociation. 
Gluon fusion arguments indicate that these states should 
be preferentially produced over negative charge conjuga- 
tion states in small- .t heavy ion collisions. 

It is of interest to speculate on the high energy limit of 
these quark model cross sections. It is apparent from the 
figures that the cross sections we find for channels that 
open at higher energies decrease in scale as the channel 
threshold increases. This is expected since the higher 
channels have a larger momentum mismatch with the 
initial state. Thus, exchanged quarks must probe the 
higher momentum region of the initial hadronic wave- 
functions. This implies that the peaks of high mass chan- 
nels will be approximately exponentially suppressed in 
Ecm, due to wavefunction suppression of the amplitudes. 
We note that final state hadronic wavefunctions do not 
affect this argument, since they are effectively averaged 
over all length scales near threshold. Furthermore, un- 
certainties due to relativistic effects or higher Fock state 
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components will not change the argument as long as the 
momentum transfer probes the confinement region of the 
wavef unctions. Uncertainties in the structure of the pion 
arising from nonperturbative spontaneous chiral symme- 
try breaking effects are also unimportant here, since the 
well-understood charmonium wavcfunctions alone suffice 
to give these general cross section features. Finally, we 
note that there is an additional suppression due to the 
nodal structure of highly excited wavefunctions. 

Once the scattering energy is large enough to probe the 
Coulombic region of the hadronic wavefunction we expect 
to find a power-law suppression of the cross section peak 
rather than an exponential suppression. This power law 
is further weakened by the nodal suppression mentioned 
above. We stress that exclusive cross section peaks must 
still fall rapidly, even in this perturbative regime. 

For two-to-two scattering the behavior of inclusive 
cross sections depends on the general behavior noted 
above and the density of states, which gives the rate 
at which new channels open with increasing ^/s. Quark 
models suggest that the density of qq resonances grows 
as a power of mass, thus the total dissociation cross sec- 
tions must decrease roughly exponentially while in the 
confinement regime, and thereafter as a power in the per- 
turbative regime. At very high energies the gluonic flux 
tube may be excited, leading to an exponential increase 
in the density of states 35] since the flux tube contains in- 
finitely many degrees of freedom. However, wavefunction 
suppression will again be exponential due to severe sup- 
pression of amplitudes containing multiply-excited string 
modes with ground state string configurations. Thus, 
for the case of two-body to two-body inclusive flavor- 
exchange scattering, nonperturbative effects cannot be 
ignored, and the cross section should decrease with in- 
creasing center of mass energy. 

The constituent quark model provides a microscopic 
foundation for the exploration of hadronic interactions at 
low energy. Thus all relevant reactions may be computed 
with the addition of no new parameters. This stands in 



contrast to effective models which must introduce new 
couplings and form factors, and which suffer from confu- 
sion over the correct degrees of freedom and dynamics to 
be employed (there are no obvious symmetries to guide 
the construction of effective lagrangians in this energy 
regime). Smu rule and pQCD computations similarly 
suffer from the notoriously poor convergence properties 
of QCD (as exemplified in renormalon ambiguities) and 
from the difficulty in extracting observables from conden- 
sates. We regard the constituent quark model as the most 
reliable tool for the investigation of these issues and in 
future plan to apply it to charmonium-nuclcon scattering 
and other reactions of interest to RHIC. 
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